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Abstract The di-neutron correlation is a spatial correlation with which two valence neutrons are 
located at a similar position inside a nucleus. We discuss possible experimental probes for the di¬ 
neutron correlation. This includes the Coulomb breakup and the pair transfer reactions of neutron-rich 
nuclei, and the direct two-neutron decays of nuclei beyond the neutron drip-line. 
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1 Introduction 

One of the most important issues in many-body physics is to clarify the nature of correlations beyond 
the independent particle picture. In nuclear physics, the pairing correlation has been well recognized 
as a typical many-boy correlation[ll [^, which leads to a characteristic even-odd staggering of biding 
energy providing an extra binding for even-mass nuclei. A pairing interaction scatters nucleon pairs 
from a single-particle level below the Fermi surface to those above, and as a consequence, each single¬ 
particle level is occupied with a fractional occupation probability. 

With the pairing correlation, one may naively expect that two nucleons forming a pair are located 
at a similar position inside a nucleus. A spatial structure of two valence neutrons has in fact attracted 
much attention in the past. One of the earliest publications on this problem is by Bertsch, Broglia, 
and Riedel, who solved a shell model for ^^°Pb and showed that the two valence neutrons are strongly 
clusterized Q . Subsequently, Migdal argued that two neutrons may be bound in a nucleus even though 
they are not bound in the vacuum [J. 
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Fig. 1 The two-particle densities for ^^Li (the left panel) and for ®He (the right panel) obtained with a 
three-body model calculation with a density dependent contact pairing interaction m- These are plotted as a 
function of neutron-core distance, ri = r2 = r, and the opening angle between the valence neutrons, ^12. The 
densities are weighted with a factor sin612. 


The strong localization of two neutrons inside a nucleus has been referred to as the di-neutron 
correlation. It has been nicely demonstrated in Ref. Q that an admixture of confiurations of single¬ 
particle orbits with opposite parity is essential to create the strong di-neutron correlation. This implies 
that the pairing correlation acting only on single-particle orbits with the same parity is not sufficient in 
order to develop the di-neutron correlation, and the pairing model space needs to be taken sufficiently 
largely so that both positive parity and negative parity states are included. 

Although the di-neutron correlation exists even in stable nuclei, it is therefore more enhanced in 
weakly bound nuclei because the admixtures of single-particle orbits with different parities are easier 
due to the couplings to the continuum spectra[l;[ 3 - Three-body model calculations have revealed that 
a strong di-neutron correlation indeed exists in weakly-bound Borromean nuclei, such as ^^Li and ®He 
iina ; [ll|; [H; [l^. For instance. Fig. [1] shows the two-particle density for the ^^Li and ®He nuclei 
obtained with the three-body model calculation with a density dependent contact pairing interaction 
m- One can see that the densities are concentrated in the region with small opening angles, that is 
nothing but the di-neutron correlation. It has been shown that the di-neutron correlation exists also in 
heavier neutron-rich nuclei [l^ : [Tsjl as well as in infinite neutron matter [T^ . The di-proton correlation, 
which is a counter part of the di-neutron correlation, has also been shown to exist in the proton-rich 
Borromean nucleus, ^'^Ne[l 7 l|. 

From these studies, the di-neutron correlation seems to have been theoretically established. How¬ 
ever, it is not straightforward to probe it experimentally. In this contribution, we discuss how one 
can probe the di-neutron correlation. To be more specific, we shall discuss the Coulomb breakup, the 
two-neutron transfer reactions, and the two-nucleon radioactivity as possible probes for the correlation. 


2 Coulomb breakup of Borromean nuclei 


Let us first discuss the Coulomb breakup reactions of Borromean nuclei, ^^Li and ®He, for which the 
experimental data have been available in Refs, [l^ fl^. Those experimental breakup cross sections, 
especially those for the ^^Li nucleus, show a strong concentration in the low excitation region, reflecting 
the halo structure of these nuclei. Moreover, the experimental data for ^^Li are consistent only with the 
theoretical calculation which takes into account the interaction between the valence neutrons, strongly 
suggesting the existence of the di-neutron correlation in this nucleus (see also Ref. 0). 

For the Coulomb breakup of Borromean nuclei, one can go one step further, given that the Coulomb 
breakup process takes place predominantly by the dipole excitation. The Coulomb breakup cross 
sections with the absorption of dipole photons are given by 


dE-^ 9 hc dE-^ 

where is the number of virtual photons, and 


dH(El) 

dE.y 


^^\{i;f\\DU,)f6{Ef-E,-E-,), 


( 1 ) 

( 2 ) 
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Fig. 2 


The geometry of a 2n halo nucleus consisting of a core nucleus and two valence neutrons. 


is the reduced E\ transition probability. In this equation, ipi and ijjj are the wave functions for the 
initial and the final states, respectively, li is the spin of the initial state, and is the operator for 
the £11 transition. For the Borromean nuclei, assuming a three-body structure with an inert core, the 
£11 operator reads [2lj |. 

^ (riYi^iri) + r2Yi^(f2)), (3) 

where the El effective charge is given by 


cei 


2Za 


(4) 


with Ac and Zc being the mass and charge numbers for the core nucleus. ri and r 2 are the coordinates 
for the valence neutrons. 

Using Eq. ([2]) and the closure relation for the final state, it is easy to derive that the total El 
strength (that is, the non-energy weighted sum rule) is proportional to the expectation value of the 
center of mass coordinate for the two valence neutrons, R^, with respect to the ground state, that is, 

i?(El)= = (5) 


with R = (ri -I- r2)l2. Even though the i3(El) strength distribution inevitably reflects both the 
correlation in the ground state and that in the final state, it is remarkable that one can extract the 
information which reflects solely the ground state properties after summing all the strength distribution. 
This implies that the average value of the opening angle between the valence neutrons can be directly 
extracted from the measured total i?(El) value once the root-mean-square distance between the valence 
neutrons, (r^„), is available (see Fig. [5]). 

This quantity is related to the matter radius and (i?^) in the three-body model as [§: . 

^ ^ ^ (rL), (6) 

where A = Ac + 2 is the mass number of the whole nucleus. The matter radii (r^) can be estimated 
from interaction cross sections. Employing the Glauber theory in the optical limit, Tanihata et al. 
have obtained = 1-57 ± 0.04, 2.48 ± 0.03, 2.32 ± 0.02, and 3.12 ± 0.16 fm for ^He, ^He, 

®Li, and ^^Li, respectively Using these values, we obtain the rms neutron-neutron distance of 

\/ {r'^n) = 3.75 ± 0.93 and 5.50 ± 2.24 fm for ®He and ^^Li, respectively. Combining these values with 
the rms core - di-neutron distance, we obtain the mean opening angle of (0„„) = 51.56l|^2 4 

and 56.2t2i’3 degrees for ®He and ^^Li, respectively [l^. These values are comparable to the result of 
the three-body model calculation, (0„„)=66.33 and 65.29 degree for ®He and ^^Li, respectively flTj . 
although the experimental values are somewhat smaller. 

An alternative way to extract the value for yj has been proposed which uses the three-body 
correlation study in the dissociation of two neutrons in halo nuclei |26l |. The two neutron correlation 
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function provides the experimental values for to be 5.9 ± 1.2 and 6.6 ± 1.5 fm for ®He, ^^Li, 

respectively [1^. Bertulani and Hussein used these values to estimate the mean opening angles and 
obtained (0nra)=83 and 66 degrees for ®He and ^^Li, respectively [^. After correcting the 
effect of Pauli forbidden transitions using the method presented in Ref. , these values are slightly 
altered to be {0„„)=74.5 and 65.2 degrees for ®He and ^^Li, respectively [^. Notice that 
these values are in a better agreement with the results of the three-body calculation [^, especially 
for the ®He nucleus. 

In the absence of the correlations, the mean opening angle is exactly (0nra)=9O degrees. The ex¬ 
tracted values of {6nn) are significantly smaller than this value both for ^^Li and ®He, providing a 
direct proof of the existence of the di-neutron correlation in these nuclei. A small drawback is that this 
method provides only the average value of and a detailed distribution is inaccessible. In reality, 

the mean opening angle is most probably an average of a smaller and a larger correlation angles in the 
density distribution, as has been shown in Fig. [TJ 


3 Two-neutron transfer reactions 

It has been recognized for a long time that two-neutron transfer reactions are sensitive to the pair¬ 
ing correlation [30j|. The probability for the two-neutron transfer process is enhanced as com¬ 
pared to a naive expectation of sequential transfer process, that is, the square of one-neutron transfer 
probability [^ . The enhancement of pair transfer probability has been attributed to the pairing 
effect, such as the surface localization of a Cooper pair [3; |^. The pair transfer reaction is thus consid¬ 
ered to provide a promising way to probe the di-neutron correlation. However, the reaction dynamics 
is rather complicated and has not yet been well established. For instance, it is only with a recent 
calculation that a theoretical calculation achieves a satisfactory agreement with the experimental data 
[^ . It would therefore be not surprising that the role of di-neutron in the pair transfer reaction has 
not yet been fully clarified. 

One such example is a relative importance of the one-step (the direct pair transfer) process and 
the two-step (the sequential pair transfer) process. In heavy-ion pair transfer reactions of stable nuclei, 
both processes are known to play a role |35| . For weakly-bound nuclei, most of the intermediate states 
for the two-step process are likely in the continuum spectra. It is still an open question how this fact, 
together with the Q-value matching condition, alters the dynamics of the pair transfer reaction of 
neutron-rich nuclei |36l | . 

On the other hand, the cross sections for the pair transfer reaction of the Borromean nuclei, ^^Li 
and ®He, have been measured recently [s^; S; |4l| . The data for the ^H(^^Li,®Li)^H reaction 

at 3 MeV/nucleon indicate that the cross sections are indeed sensitive to the pair correlation in the 
ground state of ^^Li [^. That is, the experimental cross sections can be accounted for only when 
the s-wave component is mixed in the ground state of ^^Li by 30-50%. Another important finding in 
this measurement is that significant cross sections were observed for the pair transfer process to the 
first excited state of ®Li [^ . which has made a good support for the idea of phonon mediated pairing 
mechanism [4^ . 

Further theoretical studies are apparently necessary in order to understand the connection between 
the two-neutron transfer process and the di-neutron correlation in Borromean nuclei. This has been 
left for future investigations. 


4 Two-proton and two-neutron decays of nuclei beyond the drip lines 

In the Coulomb breakup process discussed in Sec. 2, the ground state wave function of a two-neutron 
halo nucleus is firstly perturbed by the external electromagnetic field of the target nucleus. It may thus 
not be easy to disentangle the di-neutron correlation in the ground state from that in the excited states. 
The t wo-p roton radioactivity, that is, a spontaneous emission of two valence protons, of proton-rich 
nuclei [4^ is expected to provide a good tool to probe the di-proton correlation in the initial wave 
function. An attractive feature of this phenomenon is that the two valence protons are emitted directly 
from the ground state even without any external perturbation. 

Ver y re cently, the ground state two-neutron emissions have also been observed, e.g., in ^*^He [2 
HI; |4^:14^: |4§| . ^®Be ^^Li [H; [5^, and [U [H; [H; [l^. This is an analogous process of the 
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Fig. 3 The decay energy spectrum for the two-neutron emission decay of obtained with the ^“^O+n-l-n 
three-body model. For the presentation purpose, a width of 0.21 MeV has been introduced to the spectrum. 
The dotted line is for the uncorrelated spectrum, while the solid line shows the correlated spectrum for the 0^ 
state. The dashed line shows the spectrum obtained by including the pairing correlation only for the [d^/ 2 ^ 
configuration. The experimental data, normalized to the unit area, are taken from Ref. |511| . 


two-proton radioactivity, corresponding to a penetration of two neutrons over a centrifugal barrier. 
Since the long range Coulomb interaction is absent, one may hope that the ground state correlation 
can be better probed by studying the energy and the angular correlations of the emitted neutrons, as 
compared to the two-proton decays. 

FigureOshows the calculated decay energy spectrum of obtained with the '^'^0-\-n+n three-body 
model for [jj5a|. The calculations are carried out using the Green’s function method as explained 
in Refs. [^: l55l: l5fij| . together with a density-dependent contact neutron-neutron interaction, v. In 

this formalism, the decay energy spectrum is given by, 

rIP 1 

— = ^ S{E - E,) = (7) 

k 

where Ek is a solution of the three-body model Hamiltonian with energy Ek and ^ref is the wave 
function for a reference state. The reference state can be taken rather arbitrarily as long as it has an 
appreciable overlap with the resonance states of interest. Here, we employ the uncorrelated two-neutron 
state in for it with the | [ 1 ^ 3/2 0 1 ^ 3 / 2 ]^^“°^) configuration, which is the dominant configuration in 
the initial state of the proton knockout reaction of ^^F to produce In Eq. d?]), G is the correlated 
two-particle Green’s function calculated as 

G{E) = (1 + Go{E)v)-^Go{E) = Gq^E) - Go{E)v{l + Go{E)v)-^GoiE), ( 8 ) 


with the uncorrelated Green’s function, Gq, given by 


Go{E) 



€i + €2 — E - irj' 


(9) 


where 77 is an infinitesimal number and the sum includes all independent two-particle states including 
both the bound and the continuum single-particle states. To this end, we use the Woods-Saxon potential 
for the neutron -^^0 potential which reproduces the experimental single-particle energies of £ 231/2 
—4.09(13) MeV and £ 1 ^ 3/2 = 770t^g keV for [13|. The parameters for the density-dependent 
zero-range pairing interaction are determined so as to yield the decay energy of 30 keV. 

In the figure, we show the spectrum for the uncorrelated case by the dotted line. In this case, the 
spectrum has a peak at E = 1.54 MeV, that is twice the single-particle resonance energy, 0.77 MeV. 
With the pairing interaction between the valence neutrons, the peak energy is shifted towards lower 
energies. The decay spectrum obtained by including only the {d^/ 2 ]'^ configurations is shown by the 
dashed line. In this case, the peak is shifted by ^^0.5 MeV from the unperturbed peak at 1.54 MeV. 
The peak is further shifted downwards by the configuration mixing, and gets closer to the threshold 
energy. This comparison implies that the pairing correlation for the single configuration alone is not 
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Fig. 4 The angular correlation for the two emitted neutrons from the ground state decay of ^®0, that is, the 
probability distribution for the opening angle of the momentum vectors of the emitted neutrons. The solid and 
the dot-dashed lines denote the correlated and uncorrelated results, respectively. 


enough to reproduce the empirical decay spectrum, and the di-neutron correlations between the two 
neutrons also play an essential role. 

The angular distribution of the emitted neutrons can also be calculated with the two-particle 
Green’s function method [2l| . The amplitude for emitting two neutrons with spin components of 

Si and S 2 and momenta ki and k 2 reads 

fs,sAki,k 2 ) = ( 10 ) 

3,1 

where yjim is the spin-spherical harmonics, Xs is the spin wave function, and 5 is the nuclear phase 
shift. Here, M is a decay amplitude calculated to a specific two-particle final state [2l|, 

+ vGovG^ = ((jj)(°°>|(l + vGo)-^\W^), (11) 

where the unperturbed Green’s function, Gq, is evaluated at ill = ei -|- 62 . The angular distribution is 
then obtained as 

P(6»i2)=47r^ /dfcidfc2 |/sis2(fci,fci = 0,fc2,fc2 = 6 'i2)P, (12) 

Sl,S2 

where we have set z-axis to be parallel to fei and evaluated the angular distribution as a function of 
the opening angle, 0 i 2 , of the two emitted neutrons. 

Fig. U shows the angular correlation so obtained. The dot-dashed line shows the distribution ob¬ 
tained without including the nn interaction, which is symmetric around 9 i 2 = t: 12. In the presence of 
the nn interaction, the angular distribution turns to be highly asymmetric, in which the emission of 
two neutrons in the opposite direction (that is, 612 = tt) is enhanced, as is shown by the solid line. 
Grigorenko et al. have also obtained a similar result (58l |. 

This behavior reflects properties of the resonance wave function of That is, because of the 
continuum couplings, several configurations with opposite parity states mix coherently. Symbolically, 
let us write a two-particle wave function as, 

<F(r,r') = a%eir,r') + l3Wooir,r'), (13) 

where ^ee and tf'oo are two-particle wave functions with even and odd angular momentum states, 
respectively. The coefficients a and (3 are such that the interference term in the two-particle density, 
a*/3!F*g<Foo + C.C., is positive for r' = r while it is negative for r' = —r so that the two-particle density 
is enhanced for the nearside configuration with r ~ r' as compared to the far side configuration with 
r ~ —r'. This correlation appears in the opposite way in the momentum space. In the Fourier transform 
of tF(r, r'), 

^(fe,fc')= J drdr'e*^ ’’e*^'-^'tF(r,r'), 


(14) 
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there is a factor z* in the multipole decomposition of Since is +1 for even values of I and 

— 1 for odd values of I, this leads to 0;[ia 

${k,k') = aWee{k,k') - l3$oo{k,k'). (15) 

for the two particle wave function given by Eq. (11311 . If one constructs a two-particle density in the 
momentum space with this wave function, the interference term therefore acts in the opposite way to 
that in the coordinate space. That is, the two-particle density in the momentum space is hindered for 
k ^ k' , while it is enhanced for k ^ —k'. 

From this argument, we can therefore conclude that, if an enhancement in the region of 0 ^ tt in the 
angular distribution was observed experimentally, that would make a clear evidence for the di-neutron 
correlation in this nucleus, although such measurement will be experimentally challenging 

Incidentally, the tunneling decay of two fermionic ultracold atoms have been measured very recently 
[13 (see also Ref. [6l| for an application of the Gamow shell model to this phenomenon). An attractive 
feature of this experiment is that several parameters are experimentally controlable, which include the 
sign and the strength of the interaction between the particles and the shape of a decaying potential. 
It may be useful to carry out in future detailed analyses of the tunneling decay of ultracold atoms in 
order to shed light on the two-proton and two-neutron decay problems in nuclear physics. 


5 Summary 

We have discussed possible experimental probes for the di-neutron correlation in neutron-rich nuclei, 
with which two valence nucleons are located at a similar position in the coordinate space. In particular, 
we have discussed the Coulomb dissociation of Borromean nuclei, the two-neutron transfer reactions, 
and the direct two-neutron decay of the unbound nucleus. 

For the Coulomb dissociation of Borromean nuclei, even though the detailed distribution is difficult 
to extract, one can use the cluster sum rule (that is, the non-energy weighted sum rule) to deduce the 
mean value of the opening angle between the valence neutrons. We have demonstrated that the mean 
opening angle is (6(„„)=74.5 and 65.2 1 i3;q degrees for ®He and ^^Li, respectively. These values 
are signihcantly smaller than the value for the uncorrelated distribution, that is, {0nn)=9O degrees, 
clearly indicating the existence of the di-neutron correlation in these Borromean nuclei. 

For the two-neutron transfer reactions, the sensitivity of transfer cross sections to the pairing 
correlation has been well recognized for a long time. However, the reaction dynamics is rather complex, 
and the relation to the di-neutron correlation has not yet been clarified completely. Given the new 
experimental data for the two-neutron transfer reactions of the Borromean nuclei, it is now at a good 
time to discuss these reactions and clarify the reaction dynamics with the di-neutron correlation. 

For the direct two-neutron decay, we have discussed the recent experimental data of the decay 
energy spectrum for the unbound nucleus. We have shown that the decay energy spectrum can 
be accounted for only with the di-neutron correlation caused by a mixing of many configurations 
including the continuum. We have also discussed the angular correlations of the emitted two neutrons. 
We have argued that the di-neutron correlation enhances an emission of the two neutrons in the opposite 
direction (that is, the back-to-back emission), and indeed our three-body model calculation has revealed 
such feature. If the enhancement of the back-to-back emission will be observed experimentally, it will 
thus provide a direct evidence for the di-neutron correlation. 

Even though we did not discuss them in this paper, there are other possible probes for the di¬ 
neutron correlation. Those include the nuclear breakup reaction 1621 1 , th e ( p , d ) scattering at backward 
angles [H; , and the knockout reactions of Borromean nuclei jH; jH; j^; @ . It would be extremely 

intriguing if a clear and direct evidence for the di-neutron correlation could be experimentally obtained 
in near future using also these probes. 
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